State-insensitive trapping of Rb atoms: linearly versus circularly polarized lights 
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We study the cancellation of difTerential ac Stark shifts in the 5s and 5p states of rubidium 
atom using the linearly and circularly polarized lights by calculating their dynamic polarizabilities. 
Matrix elements were calculated using a relativistic coupled-cluster method at the single, double 
and important valence triple excitations approximation including all possible non-linear correlation 
terms. Some of the important matrix elements were further optimized using the experimental 
results available for the lifetimes and static polarizabilities of atomic states. "Magic wavelengths" 
are determined from the differential Stark shifts and results for the linearly polarized light are 
compared with the previously available results. Possible scope of facilitating state- insensitive optical 
trapping schemes using the magic wavelengths for circularly polarized light are discussed. Using the 
optimized matrix elements, the lifetimes of the Ad and 6s states of this atom are ameliorated. 

PACS numbers: 32.60.-|-i, 37.10.Jk, 32.10.Dk, 32.70.Cs 



I. INTRODUCTION 

Investigating the properties of rubidium (Rb) atom is 
of immense interest for a number of applications . 
It is one of the most widely used atom in quantum com- 
putational schemes using Rydberg atoms, where the hy- 
perfine states of the ground state of Rb atom are de- 
fined as the qubits It is also used to study quan- 
tum phase transitions of mixed-species with degenerate 
quantum gases [(|]. There are several proposals to carry 
out precision studies in this atom such as constructing 
ultra-precise atomic clocks 043 > probing parity non- 
conservation effects finding its permanent electric 
dipole moment [TTj etc. Also, a number of measurements 
and calculations of lifetimes for many low-lying states in 
Rb have been performed over the past few decades [T2I - 
ITsj . It is found that there are inconsistencies between the 
calculated and measured values of the lifetimes of atomic 
states in this atom 18|. In this context, it is necessary 



to carry out further theoretical studies in this atom. 

Due to simple single-core electron structure of this 
atom, it is adequate to employ advanced many-body 
methods for precise calculation of its properties which 
ultimately act as benchmark tests for the experimental 
measurements [l9l - [2l| . In this paper, we determine po- 
larizabilities of the ground 5s and excited bp states and 
study the differential ac Stark shifts between these two 
states. In this process, we also analyse the reduced ma- 
trix elements and their accuracies which are further used 
to estimate precisely the lifetimes of few excited states in 
this atom. Aim of our present study is to analyse results 
of differential ac Stark shifts from which we can deduce 
the magic wavelengths (see below for definition) that are 
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of great use in state-insensitive trapping of Rb atoms. 

Manipulation of cold and ultracold Rb atoms has been 
widely done by using optical traps [221, [l^l ■ Fo^' ^ num- 
ber of applications (such as atomic clocks and quantum 
computing [l^, HI] ) , it is often desirable to optically trap 
the neutral atoms without affecting the internal energy- 
level spacing for the atoms. However in an experimental 
set up, the interaction of an atom with the externally ap- 
plied oscillating electric field of the trapping beam causes 
ac Stark shifts of the atomic levels inevitably. For any 
two internal states of an atom, the Stark shifts caused 
due to the trap light are in general different which af- 
fects the fidelity of the experiments ^6, 27]. Katori et al. 
psj proposed the idea of tuning the trapping laser to a 
magic wavelength, "Amagic" , at which the differential ac 
Stark shifts of the transition is terminated. Using this ap- 
proach the magic wavelength for the 5s^ ^Sq — 5s5p "^Pq 
transition in ^''Sr was determined with a high precision 
to be 813.42735(40) nm [29[. McKeever et al. demon- 
strated the state-insensitive trapping of Cs at Amagic ~ 
935 nm while still maintaining a strong coupling with 
the 6si/2 — 6p3/2 transition 30]. Arora et al. [3l| calcu- 
lated the magic wavelengths for the np — ns transitions 
for other alkali atoms (from Na to Cs) by calculating dy- 
namic polarizabilities using a relativistic coupled-cluster 
(RCC) method. Theoretical values for these quantities 
were calculated at wavelengths where the ac polarizabil- 
ities for two states involved in the transition cancel. The 
data in Ref. [3l| provides a wide range of magic wave- 
lengths for the alkali-metal atoms trapped in linearly po- 
larized light by evaluating electric dipole (El) matrix 
elements obtained by linearized RCC method. In this 
paper, we try to evaluate these matrix elements consid- 
ering all possible non-linear terms in the RCC method. 
In addition, we would like to optimize the matrix ele- 
ments using the precisely known experimental results of 
lifetimes and static polarizabilities for different atomic 



2 



states and re-investigate the above reported magic wave- 
lengths in the considered atom. It is also reported in Ref. 
|3ll [ that trapping Rb atoms in the linearly polarized light 
offers only a few suitable magic wavelengths for the state- 
insensitive scheme. This persuades us to look for more 
plausible cases for constructing state-insensitive traps of 
Rb atoms using the circularly polarized light. Using the 
circularly polarized light may be advantageous owing to 
the dominant role played by vector polarizabilities (which 
are absent in the linearly polarized light) in estimating 
the ac Stark shifts. Moreover, these vector polarizabil- 
ities act as "fictitious magnetic fields", turning the ac 
Stark shifts to the case analogous to the Zeeman shifts 

This paper is organized as follows, in sections and 
mil we discuss in brief the theory of dipole polarizability 
and method used for calculating them precisely. In sec- 
tion |lVl we first discuss in detail the evaluation of matrix 
elements used for precise estimation of polarizability and 
then present our magic wavelengths first for the linearly 
polarized light following which for the circularly polarized 
light. Unless stated otherwise, we use the conventional 
system of atomic units (au), in which e, me, 47reo, and 
the reduced Planck constant h have the numerical value 
1 throughout this paper. 



II. THEORY OF DIPOLE POLARIZABILITY 



The v^^ energy level of an atom placed in a static elec- 
tric field £ can be expressed using a time-independent 
perturbation theory as [33] 
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where Ei^s are the unperturbed energy levels in the ab- 
sence of electric field, k represent the intermediate states 
allowed by the dipole selection rules and V = —D ■ £ \s 
the interaction Hamiltonian with D as the electric-dipole 
operator. Since the first-order correction to the energy 
levels is zero in the present case, therefore we can ap- 
proximate the energy shift at the second-order level for 
a weak field £ and write it in terms of dipole moments p 
as 

{P*)vk{p)kv ^2 f^s 

where SE^k = - Ek°) and = {ijjy \D\ tl^k) is 

the El amplitude between and |?/)fe) states. A more 
traditional notation of the above equation is given by 
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AEy = -^ay£^, 
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where is known as the static polarizability of the w*'' 
state which is written as 



k=^v 
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If the applied field is frequency-dependent (ac field) , then 
we can still express the change in energy as AEy = 
— ^ay£'^ with function of frequency given by 
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Since av{uj) also depends on angular momentum j and 
rrij values of the given atomic state, it is customary to ex- 
press them in a different form with mj dependent factors 
and rrij independent factors. Therefore, ay{ijj) is further 
rewritten as ISal 



a'^y{u}) + A cos Ok— al{u}) 
J 

3m2 



3 cos^ e„ - 1 



where A, 9k and 9p define degree of circular polarization, 
angle between wave vector of the electric field and z- 
axis and angle between the direction of polarization and 
z-axis, respectively. Here .4 = for the linearly polar- 
ized light implying there is no vector component present 
in this case; otherwise ^ = 1 for the right-handed and 
^ = — 1 for the left-handed circularly polarized light. In 
the absence of magnetic field (or in weak magnetic field) , 
we can choose cos{9k) = cos{9p) = 1. Here rrij indepen- 
dent factors ay, a\ and are known as scalar, vector 
and tensor polarizabilities, respectively. In terms of the 
reduced matrix elements of dipole operator they are given 
by m 
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(8) 



_SEkv+U} 8Eky~ljJ 

5jv(2jy - 1) 

6(j„ + l)(2j„ + l)(2j,+3) 

ETi I ^;]i-^y"^"'^'\(^^-\\DUk)\' 



Jk 
X 



Jk 

1 1 



SEkv + SEkv - w 



(9) 



For oj = 0, the results will correspond to the static polar- 
izabilities which clearly suggests that a}, is zero for the 
static case. 
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III. METHOD OF CALCULATIONS 

To calculate wave functions in Rb atom, we first ob- 
tain Dirac-Fock (DF) wave function for the closed-shell 
configuration [4p^] which is given by |$o)- Then the DF 
wave function for atomic states with one valence config- 
uration are defined as 

l*.) = 4l$o), (10) 

where a J, represents addition of the valence orbital, de- 
noted by V, with |$o)- The exact atomic wave function 
il'^v)) for such a configuration is determined, accounting 
correlation eflFects in the RCC framework, by expressing 

M 

= e^{l + 54|*„), (11) 

which in linear form is given by 

« {1 + T + 5J|$,). (12) 

Here T and Sy operators account excitations of the elec- 
trons from the core orbitals alone and valence orbital 
together with core orbitals, respectively. In the present 
paper, we consider Eq. pT|) instead of Ecl JIT^ as was 
taken before in our previous calculations [31|. We con- 
sider here only single, double (CCSD method) and im- 
portant triple excitations (known as CCSD(T) method 
from |$o) and !$„). 

The excitation amplitudes for the T operators are de- 
termined by solving 

($S|{i?r^}|<i>o) = 0, (13) 

where |<i>g) represents singly and doubly excited config- 
urations from |$o)- Similarly, the excitation amplitudes 
for the Sy operators are determined by solving 

{^:\{H^T}{1 + Sy}\'^>y) = {K\SMAE!i", (14) 

taking |<i>*) as the singly and doubly excited configura- 
tions from \^v)- The above equation is solved simulta- 
neously with the calculation of attachment energy Ai?°** 
for the valence electron v using the expression 

AEf* = ($,|{i?>}{l + Sy}\^y). (15) 

The triples effect are incorporated through the calcula- 
tion of Ai?°** by including valence triple excitation am- 
plitudes perturbatively (e.g. see [3^ for the detailed dis- 
cussion). 

To determine polarizabilities, we divide various corre- 
lation contributions to it into three parts as 

= atic)+a^ivc)+ai{v), (16) 

where A = 0, 1 and 2 represents scalar, vector and tensor 
polarizabilities, respectively, and the notations c, vc and 
V in the parentheses correspond to core, core- valence and 



valence correlations, respectively. The core contributions 
to vector and tensor polarizabilities are zero. 

We determine the valence correlation contributions to 
the polarizability in the sum-over-states approach [s^ 
by evaluating their matrix elements by our CCSD(T) 
method and using the experimental energies (39l-l4l| for 
the important intermediate states. Contributions from 
the higher excited states and continuum are accounted 
from the following expression 

- Cx{^y\D\¥;^), (17) 

where Cx are the corresponding angular factors for dif- 
ferent values of A and l^'i^'') is treated as the first order 
wave function to due to the dipole operator D [42| 
at the third order many-body perturbation (MBPT(3) 
method) level and given as a^(tail). Also, contributions 
from the core and core- valence correlations are estimated 
using this procedure. 

We calculate the reduced matrix elements of D be- 
tween states 1 4"/) and to be used in the sum-over- 
states approach, from the following RCC expression 

(<^>f\\{i + sl}D{i + s^}\\<^>i) 

where 'D = e^^ De^ and Ny = {^y\e'^\'^ + 

S^e"^^ Sy\^y) involve two non-truncating series in the 
above expression. Calculation procedures of these ex- 
pressions are discussed in detail elsewhere [i^, 113] . 



IV. RESULTS AND DISCUSSION 

Our aim is to determine the magic wavelengths for 
the linearly and circularly polarized electric fields for 
the 5s — 5pi/2,3/2 transitions in Rb atom. To determine 
these wavelengths precisely, we need accurate values of 
polarizabilities which depend upon the excitation ener- 
gies and the El matrix elements between the intermedi- 
ate states of the corresponding states. In this respect, 
we first present below the El matrix elements between 
different transitions and discuss their accuracies. Then 
we overview the current status of the polarizabilities re- 
ported in literature and compare our results with them. 
These results are further used to determine the magic 
wavelengths for both the linearly and circularly polar- 
ized lights. 



A. Matrix elements 

The matrix elements of Rb atom have been reported 
several times previously [ul [ll S IM!, We 
present these results from our calculations in Table H] us- 
ing the DF and CCSD(T) methods; the differences in the 
results imply the amount of correlation effects involved to 
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TABLE I: Absolute values of El matrix elements in Rb atom 
in eao using the Dirac-Fock (DP) and CCSD(T) methods. 
Uncertainties in the CCSD(T) results are given in the paren- 



theses. 


Transition 


DP 


CCSD(T) 


5si/2 ^ 


> 5pi/2 


4.819 


4.26(3) 


5Sl/2 ^ 


> 6pi/2 


0.382 


0.342(2) 


5Si/2 — 


> 7pi/2 


0.142 


0.118(1) 


5^1/2 


> 8pi/2 


0.078 


0.061(5) 


5Si/2 — 


> 9pi/2 


0.052 


0.046(3) 


5Si/2 — 


> 5p3/2 


6.802 


6.02(5) 


5Sl/2 — 


> 6p3/2 


0.605 


0.553(3) 


5Sl/2 ~ 


> 7pa/2 


0.237 


0.207(2) 


5si/2 — 


> 8p3/2 


0.135 


0.114(2) 


5Sl/2 ~ 


> 9p3/2 


0.091 


0.074(2) 


5Pl/2 ~ 


> 6Si/2 


4.256 


4.144(3) 


5Pl/2 ~ 


> 7Si /2 


0.981 


0.962(4) 


5Pl/2 ~ 


> SSi /2 


0.514 


0.507(3) 


5pi/2 - 


> 9Si/2 


0.337 


0.333(1) 


5pi/2 - 


> lOSi/2 

J./ ^ 


0.239 


0.235(1) 


5pi/2 - 


> 4^3/2 


9.046 


8.07(2) 


5Pl/2 ~ 


> 5^3/2 


0.244 


1.184(3) 


5Pl/2 ~ 


> 6^3/2 


0.512 


1.002(3) 


5pi/2 - 


> 7(^3/2 


0.447 


0.75(2) 


5pi/2 - 


> 8^3/2 


0.366 


0.58(2) 


5Pl/2 ~ 


■5> 9^3/2 


0.304 


0.45(1) 


5P3/2 ~ 


> 6Si /2 

1/ z 


6.186 


6.048(5) 


5p3/2 - 


> 7Si/2 

J./ ^ 


1.392 


1.363(4) 


5p3/2 ~ 


> 8Si/2 


0.726 


0.714(3) 


5P3/2 ~ 


> 9Si/2 


0.476 


0.468(2) 


5p3/2 ~ 


> lOSi/2 


0.338 


0.330(2) 


5P3/2 - 


> 4^3/2 


4.082 


3.65(2) 


5P3/2 - 


> 5^3/2 


0.157 


0.59(2) 


5P3/2 - 


> 6^3/2 


0.255 


0.48(2) 


5P3/2 - 


> 7^3/2 


0.217 


0.355(4) 


5P3/2 - 


> 8^3/2 


0.176 


0.272(3) 


5P3/2 - 


> 9^3/2 


0.145 


0.212(2) 


5P3/2 - 


> 4^5/2 


12.24 


10.96(4) 


5P3/2 - 


> 5^5/2 


0.493 


1.76(3) 


5p3/2 - 


> 6^5/2 


0.778 


1.42(3) 


5P3/2 - 


> 7^5/2 


0.658 


1.06(2) 


5P3/2 - 


> 8^5/2 


0.530 


0.81(1) 


5P3/2 - 


■> 9^5/2 


0.417 


0.593(5) 



evaluate these matrix elements. We also give uncertain- 
ties in the CCSD(T) results mentioned in the parenthe- 
ses in the same table. The contributions to these uncer- 
tainties come from the neglected triple excitations in the 
RCC method and from the incompleteness of the used ba- 
sis functions. The uncertainty contribution from the for- 
mer is estimated from the differences between the CCSD 
and CCSD(T) results. Some of the important matrix 
elements are determined more precisely below from the 
available experimental lifetime results of atomic states 
involving only one (strong) transition channel. However 
in case there are more than one (strong) decay channels 
associated with an atomic state, it would be intricate to 
obtain the matrix elements precisely but have been done 
by optimizing these values to reproduce the experimen- 



TABLE II: Scalar polarizability of the 5s state in Rb (in an). 
Uncertainties in the results are given in the parentheses. 



Contribution 


El amplitude 


Contribution to a.,, 


"5«l/2 








OSl/2 


~^ 5pi/2 






5Si/2 


-^■ 6pi/2 




0.361 


5Si/2 


-> 7pi/2 


0.118(1) 


0.037 


5Si/2 


^ 8pi/2 


0.061(5) 


0.009 


5si /o 


— >■ 9di It 

' ^/-^l/2 


0.046(3) 


0.005 


OSl/2 


~^ 5p3/2 


5.977(9) 


203.92(4) 


5Si/2 


6p3/2 




n Q4n 


5Si/2 


7p3/2 


0.207(2) 


0.112 


5Si/2 


^ 8p3/2 


n 1 1 /I /'o^ 
U.tt4^z ) 




5Si/2 


9p3/2 


0.074(2) 


0.013 




(c) 




9.1(5) 




(uc) 




-0.26(2) 




(tail) 




0.11(1) 


Total 






318.3(6) 



TABLE III: Dynamic polarizability of the 5s state in Rb (in 
au) at A = 1064 nm. Uncertainties in the results are given in 
the parentheses. 



Contribution 


El amplitude 


Contribution to 


"5^1/2 (^) 






5Si/2 5pi/2 


4.227(6) 


235.24(3) 


5Si/2 -> 6pi/2 


0.342(2) 


0.428 


5Si/2 7pi/2 


0.118(1) 


0.041 


5Si/2 8pi/2 


0.061(5) 


0.010 


5Si/2 9pi/2 


0.046(3) 


0.006 



5Si/2 5p3/2 


5.977(9) 


441.14(8) 


5Si/2 6p3/2 


0.553(3) 


1.114 


5Si/2 7p3/2 


0.207(2) 


0.127 


5Si/2 8p3/2 


0.114(2) 


0.035 


5Si/2 9p3/2 


0.074(2) 


0.014 


"5si/2(c) 




9.3(5) 


"5si/2(l'c) 




-0.26(2) 


"5si/2(tail) 




0.12(1) 


Total 




687.3(5) 



Experiment [52] 769(61) 



tal lifetimes in conjunction with the experimental static 
polarizabilities of different atomic states as discussed be- 
low. 

In order to evaluate the magnitude of the 5s — > 5pi/2 
El transition matrix element, we use the measured life- 
time of the 5pi/2 state which was reported as 27.75(8) ns 
in Ref. [i^. Using the fact that the 5pi/2 state decays 
only to the 5s state, the line strength of the 5s 5pif2 
transition can be obtained by combining this measured 
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TABLE IV: Scalar polarizability of the 5pi/2 state in Rb (in 
au). Uncertainties in the results are given in the parentheses. 



Contribution 


El amplitude 


Contribution to a 


"5Pl/2 \'^J 






J*l/2 ^ JPl/2 


4 997/' fi"! 




5Pl/2 ^ 6Si/2 


A 1 /''^'l 




'jPl/2 ^ ' •=1/2 


n 0iPi')(A\ 




OPl/2 ^ 0'=l/2 


n '^n7('X\ 

u.uu / J 


1 1 4 




n 'H'H'Hi^ 1 


n 4^9 


fS/ji /o - — >■ in.'ii /I 


0.235(1) 


0.215 


5pi/2 -> 4d3/2 


8.069(2) 


702.89(3) 


Kin , ^ f-\<-7 

OPl/2 — > 0"3/2 




/ .olo^^l ) 


5pi/2 6^3/2 


1.002(3) 


4.560 


5pi/2 7d3/2 


0.75(2) 


2.325(1) 


5pi/2 8^3/2 


0.58(2) 


1.320(1) 


5pi/2 9^3/2 


0.45(1) 


0.770 






9.1(5) 


a5pi/2 ivc) 




- 0.0 


Q5pi/2(tail) 




12.6(1.0) 


Total 




810.5(1.1) 



lifetime with the experimental wavelength (A = 7949.8 
A) of the corresponding transition. The value of the El 
matrix element of the 5s — > 5pi/2 transition is obtained 
from this result as 4.227(6) au. Similarly, it is possible 
to deduce the magnitude of the El matrix element of 
the 5s — > 5p3/2 transition by combining the measured 
lifetime of the 5p3/2 state, reported as 26.25(8) ns j49l |. 
with its experimental wavelength (7802.4 A). However, 
the 5p3/2 state has non-zero transition probabilities to 
the 5s and 5pi/2 states via the allowed El and the for- 
bidden Ml and E2 channels. We found from our calcula- 
tions that the transition probabilities through the forbid- 
den channels are very small and negligibly influence the 
lifetime of the 5p3/2 state; in fact lies within the reported 
experimental error bar. Neglecting these contributions, 
we extract the El matrix element of the 5s — 5p3/2 tran- 
sition to be 5.977(9) au. 

The estimated El matrix elements for the 5s — 5p tran- 
sitions from the experimental data are in close agreement 
with our calculated results within the predicted uncer- 
tainties. These results are further used, along with other 
matrix elements obtained from the CCSD(T) method, to 
calculate of polarizabilities of the 5s and 5p states. In 
Table |lll we list the the polarizability of the 5s state as 
318.3(6) au along with the detailed breakdown of the var- 
ious contributions. The most precise experimental result 
reported for this quantity as 318.79(1.42) au [50| is in ex- 
cellent agreement with our result. As shown in Table HIl 
the dominant contributions to the 5s state polarizability 
are from the 5s — 5p transitions following a significant 
contribution from the core correlation. We have calcu- 
lated core correlation contribution using the MBPT(3) 
method and the given uncertainty is estimated by scaling 



TABLE V: Scalar and tensor polarizabilities of the 5p3/2 state 
in Rb (in au). Uncertainties in the results are given in the 
parentheses. 



Contribution 


El amplitude 




al 




(v) 








5p3/2 


~^ 5Si/2 


5.977(9) 


-101.96(2) 


101.96(2) 


5p3/2 


^ DAi/2 


6.048(5) 


182.89(2) 


-182.89(2) 


5p3/2 


7Si/2 


1.363(4) 


5.036 


-5.036 


5P3/2 


8Si/2 


0.714(3) 


1.149 


-1.149 


5p3/2 


~^ 9Si/2 








5P3/2 


— i- 1 n« /„ 

^ lU^i/2 


0.330(2) 


0.215 


-0.215 


r 

OP3/2 


^ y:"3/2 


3.65(2) 


74.52(1) 


59.62(1) 


5p3/2 


^ 5^3/2 


0.59(2) 


0.988 


0.791 


5p3/2 


f Uti3/2 


0.48(2) 


0.531 


0.425 


577o /o 


— y Tdo In 

' ' 1 2 


0.355(4) 


0.264 


0.211 


^/^o/ 2 


8ci3/2 


0.272(3) 


0.147 


0.118 


5p3/2 


9^3/2 


0.212(2) 


0.086 


0.069 


5p3/2 


4^5/2 


10.89(1) 


663.4(5) 


-132.7(1) 


5p3/2 


5^5/2 


1.76(3) 


8.792(5) 


-1.758(1) 


5p3/2 


-> 6^5/2 


1.42(3) 


4.647(3) 


-0.929(1) 


5p3/2 


7^5/2 


1.06(2) 


2.353(1) 


—0.471 


5p3/2 


8d5/2 


0.81(1) 


1.304 


-0.261 


5p3/2 


9^5/2 


0.593(5) 


0.677 


-0.135 




(c) 




9.1(5) 


0.0 




(tic) 




- 0.0 


~ 0.0 


'^5P3/2 


(tail) 




13.40(1.5) 


-3.15(50) 


Total 






868.0(1.7) 


-165.9(5) 



the wave functions. Our result for the core contribution 
is in very good agreement with the result obtained using 
the random phase approximation (RPA) [F?!. Consis- 
tency in the estimated 5s polarizability value obtained 
using the 5s and 5s ^P3/2 matrix elements, 

which are obtained from the experimental lifetimes of 
the hp states, and experimental polarizability result sug- 
gests that both these matrix elements are very accurate. 
In order to test the accuracy of our results further we 
reproduce the dynamic polarizability of the 5s state at 
A = 1064 nm whose experimental value is reported as 
769(61) au [52|. As shown in Table Hill our result shows 
a large discrepancy with the experimental measurement. 
Even after replacing the above El matrix elements with 
the calculated CCSD(T) results, which are slightly larger 
in magnitude, the polarizability result still does not agree 
within the experimental error bar. Therefore, it would be 
instructive to perform another measurement of this dy- 
namic polarizability to assert this result. 

It seems from the above analysis that the calculated 
El matrix elements using the CCSD(T) method are rea- 
sonably accurate and can be further employed to ob- 
tain the polarizabilities of the 5p states. However, we 
can calculate the polarizabilities of the 5p states even 
more precisely if the uncertainties in the dominant con- 
tributing El matrix elements of the 6s — 5pi/2,3/2, 
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TABLE VI: Comparison of lifetimes (in ns) of three excited 
states in Rb atom from various theoretical and experimental 



studies. 


Level 


Reco 


Other theory" 


Expt. 


4d3/2 
44/2 
6Sl/2 


82.30(17) 
89.32(16) 
45.44(8) 


83.0(8) 
89.4(9) 
45.4(1) 


86(6)" 
94(6)" 
45.57(17)'= 



Refs: " [Hi, "Ref. [i2],"Ref. [3. 



4^3/2 5pi/2,3/2 and 4^5/2 ^P3/2 transitions are 
pushed down further. In order to do so, we evaluate the 
lifetime of the 6s state using our calculated matrix ele- 
ments as 4.144(3) and 6.048(5) in au of the 6s — 5pi/2 and 
6s — 5p3/2 transitions, respectively. We obtain its lifetime 
as 45.44(8) ns against the experimental result 45.57(17) 
ns [Ti] with branching ratios 34% to the 5pi/2 state and 
66% to the 5p3/2 state neglecting the observed insignif- 
icant transition probabilities to the 5s and 4d states. 
Since we are able to obtain more precise lifetime for the 
6s state using our calculated matrix elements than the 
measurement, we assume these calculated El matrix el- 
ements are more precise than what we would have ob- 
tained from the known experimental lifetime result. 

We would further like to use the above El matrix ele- 
ments to produce the experimental polarizability of the 
5pi/2 state from which we anticipate to estimate the El 
matrix element of the 5pi/2 — ^ 4^3/2 transition accu- 
rately. Since no direct measurement of the polarizability 
of the 5pi/2 state is known to us from the literature, we 
use the differential polarizability of the 5s — >■ 5pi/2 transi- 
tion which is reported as 492.20(7) au [HI]. In fact, using 
the differential polarizability here is advantageous for the 
following three reasons: (i) we have already determined 
the polarizability of the 5s state precisely, (ii) the differ- 
ential polarizability is not affected by the uncertainty of 
the core correlation contribution, and (iii) precise values 
of few important matrix elements contributing towards 
the 5p state polarizability are known to some extent from 
the above analysis. By adding the experimental differen- 
tial polarizability with the precisely known polarizability 
of the 5s state we consider the experimental 5pi/2 state 
polarizability as 810.6(6) au; indeed this result will not 
meddle the above advantages. We find from our calcu- 
lations that the El matrix elements of the 5pi/2 — 5s, 
5pi/2 — 6s, and 5pi/2 — 4c?3/2 transitions have crucial con- 
tributions to the 5pi/2 state polarizability. Substituting 
the precisely known values of the first two elements to re- 
produce the experimental 5pi/2 polarizability result, the 
El matrix element of the 5pi/2 transition is set 

as 8.069(2) au. This agrees well with our calculated result 
8.07(2) au. From this analysis, we estimate theoretical 
value of the 5pi/2 state polarizability to be 810.5(1.1) 
au; contributions from various parts are given explicitly 
in Table HVl 

It can be noticed that the 4c?5/2 state has only one al- 
lowed decay channel to the 5p3/2 state. Therefore if the 



lifetime of the state is known precisely then the El 

matrix element of the 5p3/2 ^^5/2 transition can be es- 
timated accurately from this data. There are two experi- 
mental results for the lifetime of the 4^5 /2 state reported 
as 89.5 [54] and 94(6) ns [l2]- From the former result 
which is the latest, we deduce the El matrix element of 
the above transition to be about 10.89 au which reason- 
ably agrees with our CCSD(T) result 10.94(6) au. Matrix 
element obtain from the later lifetime data gives much 
lower absolute value with very large uncertainty com- 
pared to our calculated result which is not of our interest. 
Similarly, the lifetimes of the 4(j3/2 state are reported 
as 83.4 ns ^ and 86(6) ns The 4^/2 state has 

two strong allowed decay channels to the 5pi/2 and 5p3/2 
states. By combining the above El matrix element for 
the 5pi/2 — >■ 4(^3/2 transition and the lifetime of the 4^3/2 
state as 83.4 ns (the reason for not considering the other 
value is same as cited above), we predict the El matrix el- 
ement of the 5p3/2 — >■ 4(^3/2 transition to be about 3.5 au. 
In order to find this matrix element more precisely, we use 
the experimental results for the scalar and tensor polariz- 
abilities of the 5p3/2 state which are reported as 857(10) 
and —163(3) in au [IB], respectively. Our calculation 
shows that the major contributions to these polarizabili- 
ties come from the matrix elements of the 5p3/2 5s 1/2, 
5p3/2 6si/2, 5p3/2 -J> 4^3/2 and 5^3/2 4(1^/2 transi- 
tions. From the sensitivity in the given precision of the 
El matrix element of the 5p3/2 — > 4^3/2 transition to be 
able to reproduce the scalar and tensor polarizabilities in 
their respective error bars, we get a lower bound for this 
matrix element as 3.6 au. Without any loss of quality, 
we retain our CCSD(T) result, i.e. 3.65(2) au, as the 
most precise value for this matrix element. Using these 
optimized results and the combined experimental values 
of the scalar and tensor polarizabilities of the 5p3/2 state, 
we get the best value for the El matrix element of the 
5p3/2 — ^ 4(^5/2 transition to be 10.89(1) au. 

Now we list below the optimized El matrix elements 
(in au) obtained from the above analysis apart from our 
calculated results as: 

(5s]]£>]]5pi/2) = 4.227(6) 

(5s]]D]]5p3/2) = 5.977(9) 

(5pi/2]]i?]]4d3/2) = 8.069(2) 

{5p3/2\\D\\Ad5/2) - 10.89(1). (19) 

B. Lifetimes of few excited states 

Since we are now able to estimate some of the El ma- 
trix elements more precisely than the previously known 
results, we would like to use them further to estimate 
lifetimes of first few excited states in Rb atom accu- 
rately. The matrix elements for the (5s]]r']]5pi/2) and 
(5s] ]D]] 5^)3/2) transitions were obtained from the lifetime 
measurements, so we still consider the most accurately 
known lifetimes of the 5pi/2 and 5p3/2 states from the 
experiment as 27.75(8) ns and 26.25(8) ns, respectively. 
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FIG. 1: (color online) Magic wavelengths identified by arrows 
for the 5pi /2 — 5s transition in Rb using the linearly polarized 
hgh 




600 700 800 900 1000 1100 1200 1300 1400 1500 

wavelength (nm) 

FIG. 2: (color online) Magic wavelengths identified by arrows 
for the 5p3/2 ^ 5s transition in Rb using the linearly polarized 
hght. 



We now determine the lifetimes of the 4d and 6s states 
using the El matrix elements listed in Eq. (jl9|) and from 
our calculations which are given in TableHl The estimated 
lifetimes are mentioned in Table IVTl as recommend (Reco) 
values and compared with the other available experimen- 
tal and theoretical results in the same table. 



C. Status of the polarizability results 



TABLE VII: Comparison of the static and ac polarizabilities 
(in au) in Rb atom for the 5s, 5pi/2 and 5p3/2 states with 
other experiments and theory. 













Present 


318.3(6) 


810.5(1.1) 


868.0(1.7) 
867'* 


-165.9(5) 


Other 


317.39" 


805^* 


-167^* 


Other 


318.6(6)'' 


807" 


870*= 


-171"= 


Exp. 


318.79(1.42)^= 


810.6(6)^ 


857(10)» 


-163(3)5 



References: " [62] 
' [63] 
" [50] 
[31| 

^ M 

f [53, 65] 
» [55] 



TABLE VIII: Magic wavelengths Amagic for the linearly po- 
larized light above 600 nm for the 5pi/2 — 5s and 5p3/2 — 5s 
transitions in Rb and the corresponding values of polarizabil- 
ities at the magic wavelengths. The wavelengths (in vacuum) 
are given in nm and polarizabilities are given in au. The given 
mj values correspond to the 5p states. 



Transition: 5pi 



/2 



5s 







AmagicRef. [31] 


Q^('^magic) 


1/2 


606.25(1) 


606.2(1) 


-443.3 


1/2 


618.09(2) 


617.7(7) 


-490 


1/2 


727.343(2) 


727.35(1) 


-1876 


1/2 


761.6221(2) 


761.5(1) 


-5270 


1/2 


787.633(2) 


787.6(1) 


5382 


1/2 


1350.801(9) 


1350.9(5) 


475.5 


Transition: 


5p3/2 - 5s 






I'm A 


-^magic 


AmagicRef. [31] 


'^('^magic) 


1/2 


614.70(1) 


614.7(1) 


-477 


3/2 


626.62(3) 


626.2(9) 


-529 


1/2 


627.70(1) 


627.3(5) 


-534 


1/2 


740.063(2) 


740.07(1) 


-2493 


1/2 


775.868(1) 


775.84(1) 


-20030 


3/2 


775.8228(2) 


775.77(3) 


-19917 


3/2 


790.018(2) 


789.98(2) 


53 


1/2 


792.022(1) 


792.00(1) 


-6973 


1/2 


1414.83(3) 


1414.8(5) 


455 



To affirm the broad interest of studying polarizabili- 
ties in Rb atom, we discuss briefly below about various 
experimental and theoretical works in the evaluation of 
polarizabilities of the 5s and bp states reported so far 
in the Rb atom. There were several measurements car- 
ried out on Stark shifts in Rb atom almost two decades 
ago [H, [sol - issj from which the polarizabilities of the 
5s ground state and few excited states were estimated. 
Hunter and coworkers had observed the dc Stark shifts 



of the Dl line in Rb using a pair of cavity stabilized 
diode lasers locked to resonance signals [H, [H, [s^l . In 
another work. Tanner and Wieman had used a crossed- 
beam laser spectroscopy with frequency stabilized laser 
diodes to measure the differential Stark shift of D2 line 
[ssj . Marrus et al. had used a atomic beam method 
long ago to measure the Stark shift from which both the 
scalar and tensor polarizabilities of the state were 

determined [SQi] . 
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The extensive calculation of polarizabilities in Rb atom 
was first carried out by Marinescu et al. using an l- 
dependent model potential [1^1 ■ In this work, the infinite 
second-order sums in the polarizability calculations were 
transformed into integrals over the solutions of two cou- 
pled inhomogeneous differential equations and the inte- 

f'als were carried out using Numerov integration method 
0]. In 2004, Zhu et al. employed the RCC method 
to calculate the scalar and tensor polarizabilities of the 
ground and the first p excited states in alkali atoms [13] . 
The results obtained using the RCC method were sub- 
stantially improved over the earlier calculations based on 
the non-relativistic theories. Later, Arora et al. extended 
these calculations to obtain frequency-dependent scalar 
and tensor polarizabilities of the ground and first excited 
5p states in Rb [3l|, H^, HH using the RCC method at the 
linearized singles, doubles and partial triples excitations 
level (SDpT method). 

As has been discussed earlier, we have optimize at least 
seven important El matrix elements which are crucial 
in obtaining the polarizabilities of the 5s and 5p states 
and the other matrix elements have been obtained using 
the CCSD(T) method which includes all the non-linear 
terms. Therefore, the predicted polarizabilities of the 5s 
and 5p states obtained in this work are expected to be 
accurate enough to employ them further in the determi- 
nation of the magic wavelengths in Rb atom, which is 
the prime motivation of the present work. In Table IVIH 
we compare our polarizability results with the other re- 
ported values. Our results are more polished over the 
earlier studied results mainly due to the optimization of 
the matrix elements. 

In our knowledge, there are no experimental and/or 
theoretical results on vector polarizabilities of the 5s and 
5p states for any wavelengths available in Rb atom to 
compare with the present calculations. Accuracy in these 
polarizabilities will determine correct values of the magic 
wavelengths for the circularly polarized light in this atom. 
Since the El matrix elements required to determine the 
vector polarizabilities are same as required for the calcu- 
lation of scalar and tensor polarizabilities, we expect a 
similar precision in our used vector polarizability results 
as discussed in the last subsection. We shall present the 
vector polarizabilities in the 5s and bp states at a given 
wavelength (say close to a particular Amagic value) so that 
if necessary our results can also be further verified by any 
other study. 



D. AC Stark shifts and magic wavelengths 

Following Eq. © , the ac Stark shift A£'„ of an atomic 
energy level Ey due to the external applied ac electric 
field £, in the absence of any magnetic field, can be 



(a„) polarizabilities at A=770 nm in Rb. Uncertainties in the 
results are given in the parentheses. 



Contribution 




q1 








5Si/2 — >■ 5pi/2 






5Si/2 — >■ 6pi/2 


U.OlO 


— U.ODO 


5Si/2 7pi/2 


0.047 


0.044 


5Si/2 8pi/2 


0.011 


0.010 


5Si/2 9pi/2 


0.006 


0.005 




—7615(1) 


—7716(1) 


5Si/2 -> 6p3/2 


1.339 


0.731 


5Si/2 -> 7p3/2 


0.144 


0.067 


5Si/2 -> 8p3/2 


0.039 


0.017 


5Si/2 ^> 9p3/2 


0.016 


0.007 


05^1/2 (c) 


9.2(5) 


0.0 




-0.26(2) 


~ 0.0 


Otail 


0.14(1) 


0.002(1) 


Total 


-9180(1.1) 


-4462(1.1) 



parametrized in terms of ao, o-i and a2 as 32] 

1 777 . 

2 Jv 



3m| - + 1) 
jvi2jy - 1) 



(20) 



In this formula, the frequency tu is assumed to be several 
line-widths off-resonance. The differential ac Stark shift 
for a transition is defined as the difference between the 
Stark shifts of individual levels. For instance, the inter- 
ested differential ac Stark shifts in our case are for the 
5pi — 5s transitions (with i = 1/2,3/2) which are given 
by 



AE^p, 
1 



AEi 



5s 



(21) 



where we have used the total polarizabilities of the re- 
spective states. Since the external electric field £ is arbi- 
trary, we can verify the frequencies or wavelengths where 
«5pi = Q^5s, for the null differential ac Stark shifts. 

In order to estimate the total polarizability for any 
particular set of jv and mj values, we need to deter- 
mine the scalar, vector and tensor polarizabilities. Magic 
wavelengths are calculated for a continuous values of fre- 
quencies (can also be expressed in terms of wavelength 
A) by plotting the total polarizability for different states 
against the A values. The crossing between the two po- 
larizabilities at various values of wavelengths will corre- 
spond to Amagic ■ Trapping of Rb atoms is convenient at 
these wavelengths as was stated in the beginning. As 
pointed out in [sij . the linearly polarized lattice scheme 
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TABLE X: Contributions to the 5pi/2 scalar (a^l) and vector 
(ai) polarizabilities at A=770 nm in Rb. Uncertainties in the 
results are given in the parentheses. 



Contribution 


n 

aS 


1 




(v) 






5pi/2 


5pi/2 


1567.1(2) 


3254.4(4) 


5pi/2 


-> 6Si/2 


-85.029(5) 


292.38(2) 


5jDi /o 

i-j A 


7Si/2 


46.676(4) 


-88.283(7) 


5pi/2 


^ 8Si/2 


3.020 


-4.764 


5pi/2 


9Si/2 


0.954 


-1.382 




— ^ 1 '^1 
?■ -LUAi/2 


0.412 


-0.570 


5pi/2 


-5> 4^3/2 


-262.99(1) 


—504.00(2) 


5pi/2 


5ci3/2 


382.94(2) 


379.02(2) 


5pi/2 


-5> 6^3/2 


13.029(1) 


10.504(1) 


5pi/2 


^ 7^3/2 


5.035(2) 


3.694(2) 


5pi/2 


-^> 8^3/2 


2.565(1) 


1.787(1) 


5pi/2 


9d3/2 


1.413 


0.954 




(c) 


9.2(5) 


0.0 




(vc) 


~ 0.0 


~ 0.0 


at ail 




17.6(20) 


3.8(4) 


Total 




1711(2) 


3347.7(4) 



offers only a few cases in which the magic wavelengths 
are suitable from the experimental point of view. There- 
fore, we would like to explore the idea of using the circu- 
larly polarized light for which the magic wavelengths need 
to be determined separately for each magnetic quantum 
number rrij. 

In the next two subsections, we shall discuss about 
the magic wavelengths for the 5pi/2,3/2 ~ 5s transitions 
for both the linearly and circularly polarized lights. The 
reason for bringing up the issue of magic wavelengths for 
the linearly polarized lights in these transitions is that 
since we have obtained the most accurate results for all 
the static polarizabilities it is expected that we will get 
better results for the magic wavelengths using our opti- 
mized set of El matrix elements. This will also help us in 
making a comparison study between the results obtained 
from the linearly and circularly polarized lights. 



E. Case for the linearly polarized optical traps 

Since we are interested in optical traps and the pre- 
vious study [m reveals that the magic wavelengths for 
the 5s — 5p transitions at which the Rb atom can be 
trapped using the linearly polarized lights lie in between 
600 — 1500 nm, we try to find out the null differential 
polarizabilities in this region. In Fig. [U we plot the to- 
tal polarizabilities due to the linearly polarized lights for 
both the 5s and 5pi/2 states. As seen in the figure, the 
5s state dynamic polarizabilities are generally small in 
this region except for the wavelengths in close vicinity 
to the 5s — 5pi/2 resonance (at 795 nm) and 5s — 5ps/2 
resonance (at 780 nm). However, the 5pi/2 state has sev- 



TABLE XI: Contributions to the 5p3/2 scalar (aJJ) and vector 
(al) and tensor (a^) polarizabilities at A=770 nm in Rb. 
Uncertainties in the results are given in the parentheses. 

Contribution aJJ al 



5p3/2 


5Si/2 


3807.5(7) 


11575(2) 


-3807.5(7) 


5p3/2 


-5> 6Si/2 


-85.017(9) 


452.76(5) 


85.017(9) 


5p3/2 


>. 7-„ 

— ^ 'Si/2 


dQ ■] Q-\ f (i\ 


— lyD.oOl^z ) 


ao ■] Q-\ f (i\ 


5p3/2 


— 8Si/2 


Q 1 O/l 


— / .Db {[I) 


Q 1 fiA 

— o.iy4 


5p3/2 


9Si/2 


n QQ/1 
u.yo4 


— Z.lDo 


— u.yo4 


5p3/2 


1 1 Ac. 




— U.ooO 




5p3/2 


4d3/2 


-25.311(5) 


60.32(1) 


-20.249(4) 


5p3/2 


-S> 5^3/2 


-61.57(1) 


74.47(1) 


-49.25(1) 


5p3/2 


^ 003/2 








5P3/2 


^ 7c?3/2 


0.591 


—0.527 


0.473 


5p3/2 


— V Rfii 
? OLt3/2 


0.293 


—0.248 


0.234 


5z?^ /o 


— >• 9(^3/2 


0.162 


-0.133 


0.130 


5p3/2 


^ 44/2 


-225.3(2) 


-805.4(6) 


45.06(4) 


5p3/2 


-5> 5^5/2 


-564.1(3) 


-1023.3(6) 


112.8(1) 


5p3/2 


-> 6^5/2 


14.06(1) 


20.70(1) 


-2.811(2) 


5p3/2 


-> 7^5/2 


5.264(2) 


7.046(3) 


-1.053 


5p3/2 


^ 84/2 


2.597(1) 


3.298(1) 


-0.519 


5p3/2 


-5> 94/2 


1.2700 


1.562 


-0.254 




(c) 


9.3(5) 


0.0 


0.0 




(vc) 


~ 0.0 


~ 0.0 


~ 0.0 


Cftail 




19(2) 


6.9(7) 


-4.7(9) 


Total 




2973(2) 


10163(5) 


-3714(1) 



eral resonances in the considered wavelength range. It is 
generally expected that the 5pi state polarizability will 
cross the 5s state polarizability in between each pair of 
resonances. We found total six magic wavelengths for the 
5pi/2 ~^ 5s transition in between the five resonances. 

However, the case for the 5p3/2 — >■ 5s transition is dif- 
ferent owing to the presence of non-zero tensor contri- 
bution of the 5p3/2 state. As shown in Fig. [2l we get 
different magic wavelengths for the 5p3/2 — ^ 5s transi- 
tion at TOj — ±1/2 and = ±3/2 sub-levels of the 
5p3/2 state. There are few wavelengths in between reso- 
nances where 0:5^3^2 with rrij = ±3/2 contribution is not 
same as the as^. This leads to reduction in the number 
of magic wavelengths for this transition. For example, 
we did not find any Amagic between the 5p3/2 — 4^3/2,5/2 
resonances (at 1529 nm) and the 5p3/2 —6s resonance (at 
1367 nm) for mj = ±3/2 sublevels of the 5p3/2 state. 

We have limited our search for the magic wavelengths 
where the differential polarizabilities between the 5s and 
hpj states are less than 0.5%. Based on all these data, 
we list now Amagic (in vacuum) above 600 nm in Table 
IVIIII for the 5pi/2 — 5s and 5p3/2 — 5s transitions in Rb 
atom and compare them with the previously known re- 
sults. The present results are improved slightly due to 
the optimized El matrix elements used here. The un- 
certainties in our magic wavelength results are found as 
the maximum differences between the ± Jass and 
asp ± ^Qfsp contributions with their respective magnetic 
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FIG. 3: (color online) Magic wavelengths identified by arrows 
for the 5pi/2 — 5s transition in Rb using the left-handed cir- 
cularly polarized light. 



quantum numbers, where the Sa are the uncertainties in 
the polarizabihties for their corresponding states. 

The reason for not acquiring sufficient number of magic 
wavelengths for the 5p3/2 — 5s transition lies in the fact 
that extra contribution from the tensor polarizability to 
the total 5p3/2 polarizability is not compensated by the 
counter part of the 5s state. The idea of using the cir- 
cularly polarized light to obtain magic wavelengths for 
the 5p3/2 — 5s transition is triggered from that fact that 
the extra contribution from the tensor polarizability to 
the 5p3/2 state might be cancelled by the vector polariz- 
ability contributions or the vector polarizabihties are so 
large that they may play a dominant role in determining 
the differential polarizabihties. This would be evident in 
the following subsection. 



F. Case for the circularly polarized optical traps 

As mentioned previously, polarizabihties for the circu- 
larly polarized light have extra contribution from the vec- 
tor component of the tensor product between the dipole 
operators. This extra factor is expected to provide better 
results for state-insensitive trapping. First, we present 
the scalar, vector and tensor dynamic polarizabihties of 
the 5s, 5pi/2 and 5p3/2 states in Tables HXl 1x1 and IXI[ 
respectively, at A = 770 nm to perceive their general be- 
havior. The choice of this wavelength is deliberate for 
being close to one of the magic wavelengths for the cir- 
cularly polarized light (e.g. see Table (jXIII) and (|XIIII) '). 
Hereafter we shall consider the left-handed circularly po- 
larized light for all the practical purposes as the results 
will have a similar trend with the right-handed circularly 
polarized light due to the linear dependency of degree of 
polarizability A in Eq. ([20]) . Nevertheless, the left or 
right handed polarization in the experimental set up is 
just a matter of choice. 



TABLE XII: Magic wavelengths Amagic above 600 nm for the 
5pi/2 — 5s transition in Rb and the corresponding values of 
total polarizabihties at the magic wavelengths for the left- 
circularly polarized laser beam. The wavelengths (in vacuum) 
are given in nm and polarizabihties are given in au. The given 
mj values are for the 5p states. 

Transition: 5pi/2 — 5s 



Amagic 



CK(Amagic) Amagic (avg 



1/2 600.83(14) 
-1/2 607.98(1) 



-1/2 

1/2 
-1/2 

-1/2 

-1/2 



616.77(2) 

721.628(23) 
728.843(1) 

761.176(1) 

1306.08(1) 



-405 
-428 

-461 

-1449 
-1633 

-3424 

504 



604(7) 



617 



725(7) 



761 



1306 



For the sake of completeness of our study, we also 
search for magic wavelengths in the 5s — 5pi/2 transition 
in Rb atoms using the circularly polarized light although 
a fairly large number of magic wavelengths for this tran- 
sition is found using the linearly polarized light. For this 
purpose, we plot net dynamic polarizability results of the 
5s and 5pi/2 states in Fig. [3] using the circularly polar- 
ized light against different values of wavelength. The 
figure shows that the total polarizability of the 5s state 
for any values of A is very small except for the wave- 
lengths close to the two primary resonances. Due to the 
rrij dependence of the vector polarizability coefficient in 
Eq. (PO)) . the crossing occurs at a different wavelength 
for the different values of rrij in between two 5pi /2 reso- 
nances. As shown in Table IXlIl we get set of five magic 
wavelengths in between seven 5pi/2 resonances lying in 
the wavelength range 600-1400 nm. Out of these five sets 
of magic wavelengths three sets of the magic wavelengths 
occur only for negative values of mj. Thus, the number 
of convenient magic wavelengths for the above transition 
is less than the number of magic wavelengths obtained 
for the linearly polarized light. This advocates for the 
use of linearly polarized light in this transition, though 
choice of the circularly polarized light is not bad at all. 
The rrij dependence of traps and the difficulties in build- 
ing a viable experimental set up in the case of circularly 
polarized light could be the other major concern. 

In this work, we also propose the use of "switching 
trapping scheme" (described below) which may solve the 
problem in cases where state-insensitive trapping is only 
supportive for the negative mj sublevels of 5p states. We 
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TABLE XIII: Magic wavelengths Amagic above 600 nm for 
the 5p3/2 — 5s transition in Rb and the corresponding values 
of total polarizabilities at the magic wavelengths for the left- 
circularly polarized laser beam. The wavelengths (in vacuum) 
are given in nm and polarizabilities are given in au. The given 
rrij values are for the 5p states. 



Transition: 5p: 



'3/2 



5s 





-^magic 






1/2 


613.25(3) 


-447 




-1/2 


615.51(1) 


-456 


616(5) 


-3/2 


618.15(2) 


-466 




3/2 


630.142(1) 


-516 




1/2 


628.30(1) 


-508 










628(5) 


-1/2 


626.95(1) 


-502 




-3/2 


625.04(3) 


-494 




3/2 


746.737(15) 


-2328 




1/2 


738.794(32) 


-1964 










742(8) 


-1/2 


740.587(1) 


-2037 




-3/2 


742.262(1) 


-2109 




3/2 


775.836(5) 


-6231 




1/2 


775.834(7) 


-6230 










775.8(2) 


-1/2 


775.789(3) 


-6215 




-3/2 


775.693(2) 


-6183 




1/2 


783.883(13) 


-10925 




-1/2 


787.547(4) 


-16431 


786(4) 


-3/2 


776.497(4) 


-16318 




1/2 


1454.4(9) 


453 




-1/2 


1387.1(1) 


473 


1382(149) 


-3/2 


1305.9(1) 


504 





observed that the same magic wavelength will support 
state- insensitive trapping for negative mj sublevels if we 
switch the sign of A and mj of 5s state. In other words, 
the change of sign of A and mj sublevels of 5s state will 
lead to the same result for the positive values of mj sub- 
levels of 5p states. 

Here we give more emphasis on finding more magic 
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FIG. 4: (color online) Magic wavelengths identified by arrows 
for the 5p3/2 — 5s transition in Rb using the left-handed cir- 
cularly polarized light. 



wavelengths for the 5s — 5p3/2 transition which can be 
used in the state-insensitive trapping scheme for the Rb 
atom. In Table IXIIIi we list a number of Amagic for the 
5s — 5p3/2 transition in the far-optical and near infrared 
wavelengths along with the uncertainties in the Amagic 
and the polarizabilities at the Amagic values. We also list 
the Amagic (avg) values in the table which are the average 
of the magic wavelengths at different rrij sublevels. The 
error in the Amagic (avg) is calculated as the maximum 
difference between the magic wavelengths from different 
rrij sublevels. For this transition we get a set of six magic 
wavelengths in between seven 5p3/2 resonances lying in 
the wavelength range 600-1400 nm (i.e. 5p3/2 — idj res- 
onance at 1529 nm, — 6s resonance at 1367 nm, 
5p3/2 — 5s resonance at 780 nm, 5p3/2 — 5dj resonance at 
776 nm, 5p3/2 — 7s resonance at 741 nm, 5p3/2 — 6c?j res- 
onance at 630 nm, and 5p3/2 — 8s resonance at 616 nm). 
Five out of six magic wavelengths support a blue detuned 
trap (predicted by the negative values of dynamic polar- 
izability). Out of these five magic wavelengths the magic 
wavelength at 628 nm and 742 nm are recommended for 
blue detuned traps. The magic wavelength at 742 nm 
supports a stronger trap (as shown by a larger value of 
the polarizability at this wavelength in Fig.(|3])). The 
magic wavelength at 775.8 nm is very close to the res- 
onance and might not be useful for practical purposes. 
The magic wavelength at 1382 nm supports a red de- 
tuned optical trap. It can be observed from Table IXlTTl 
that mj = 3/2 sublevel does not support state-insensitive 
trapping at this wavelength. However, using a switching 
trapping scheme as described in the previous paragraph 
will allow trapping this sublevel too. The magic wave- 
length at 1382 nm is recommended owing to the fact 
that it is not close to any atomic resonance and supports 
a red-detuned trap which was not found in the linearly 
polarized trapping scheme. 



12 



V. SUMMARY 

In conclusion, we have employed the relativistic cou- 
pled cluster method in the singles, doubles and triples ex- 
citations approximation to determine the electric dipole 
matrix elements in rubidium atom. Some of the impor- 
tant matrix elements were further optimized using the 
experimental lifetimes of few excited states and static po- 
larizabilities of the ground and 5pi/2,3/2 excited states. 
These optimized matrix elements were then used to im- 
prove the precision of the available lifetime results for 
some of the low-lying excited states in the considered 
atom. We also observe disagreement between our calcu- 
lated dynamic polarizability with a measurement at the 
wavelength 1064 nm using the above optimized matrix 
elements. 

We have compared the static and dynamic polarizabil- 
ity results from various works and reported the improved 
values of the magic wavelengths for the 5s — 5pi/2 tran- 
sition using the linearly polarized light. Issues related 
to state-insensitive trapping of rubidium atoms for the 



5s — 5p3/2 transition with linearly polarized light are 
discussed and use of the circularly polarized light is em- 
phasized. Finally, we evaluate six set of magic wave- 
lengths for the 5s — > 5p^/2 transition which can be used 
for the above purpose out of which we have recommended 
two magic wavelengths at 628 nm and 742 nm for the 
blue detuned optical traps and 1382 nm for the red de- 
tuned optical traps. We also proposed the use of a switch- 
ing trapping scheme for the magic wavelengths at which 
the state-insensitive trapping is supported only for either 
positive or negative mj sublevels of bp states. 
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